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ABSTRACT 

One of the more significant technological problems associated with 
the orbital operation of large astronomical telescopes is the fabrication 
and maintenance of the primary mirror surface to the tolerance required 
for diffraction-limited performance. An interesting approach to the 
solution of this problem involves continuously measuring and automa- 
tically correcting the optical surface of a thin deformable mirror by 
means of discrete actuators located on its rear surface. The real- 
ization of diffraction-limited performance from a telescope in space 
by this method rests on the ability of the designer to achieve extremely 
accurate control of a highly complex, interacting, multivariable system. 
This paper presents the results of a detailed study of the discrete 
control of linear distributed systems with specific application to 
the design of a practical controller for a plant representative of a 
telescope primary mirror for an orbiting astronomical observatory. 

The problem of controlling the distributed plant is treated by 

employing modal techniques to represent variations in the optical 

figure. Distortion of the mirror surface, which arises primarily from 

thermal gradients, is countered by actuators working against a backing , 

structure to apply a corrective force distribution to the controlled 

surface. Each displacement actuator is in series with a spring attached 
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to the mirror by means of a pad intentionally introduced to restrict 
the excitation of high-order modes. Control is then exerted over 
a finite number (equal to the number of actuators) of the most 
significant modes. 

Through the application of the modal expansion technique the mirror 
equation of motion is transformed to a set of uncoupled, linear, time- 
invariant, ordinary differential equations. The desired dynamic 
response and static accuracy may then be achieved by the application of 
classical single -variable design techniques. The formulation of a 
quadratic performance index which incorporates a measure of image quality 
permits determination of the trade-off between the number of actuators 
and optical purity. A criterion for defining actuator placement and 
pad size is presented which minimizes the tendency of the controller to 
excite the unmonitored modes. 
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I. INTRODUCTION 


The class of distributed plants considered in this paper is 

restricted to those described by linear, time -invariant, separable 

operators where control is derived from a finite number of discrete 

inputs. Application of the modal expansion approach^- 1 -) converts the 

distributed-parameter problem to one of a multivariable nature which 

( 2 ) 

readily yields to decoupling techniques. ' ' Classical single-variable 
control methods are applied to decoupled system dynamics defined in 
terms of the eigenvalues of the linear operator whose eigenfunctions 
are assumed to form a complete orthonormal set. While the results 
are applicable to the general problem of controlling linear distributed 
systems, the motivation for the study is a direct result of its 
relevance to one of the central problems of large orbiting telescope 
technology . 

Elimination of the effects of the earth's atmosphere give orbiting 
telescopes significant advantages relative to earth-based telescopes. 
Figure 1 displays percent transmission of the incident radiation 
through the earth's atmosphere as a function of wavelength. In 
the portion of the spectrum shown, the earth's atmosphere is opaque 
to radiation shorter than 3 X 10” ^ cm and longer than 3 X 10"^ cm, with 
numerous gaps between these extremes. Since the entire spectrum of 
radiated energy becomes available for study in an atmosphere -free 
environment the additional spectral coverage would permit studies 

. . (4) 

involving galactic nebulae and cool stars (stellar evolution). 
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Figure 1.- Transmission of incident radiation through the earth's atmosphe 
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Secondly, refraction anomalies associated with the turbulence of the 
earth's atmosphere limit the resolving power of earth-based telescopes; 
consequently, the 200- inch (508 cm) Hale telescope at Mt. Palomar has 
no better resolution than a high quality telescope of approximately 
15 inches (38 cm) diameter. ' Placing a 120-inch (3 meter) telescope 
in orbit would yield an increase in resolution of at least a factor of 
7 relative to land-based telescopes and 3 relative to present space 
telescopes. This significant improvement is useful in studying double 

(4) 

stars and in planetary photography. Figure 2 shows two views of 
the same portion of the lunar surface. The photograph on the right 
was taken by the 120-inch (3 meter) reflector at the Lick Observatory 
and represents about the best resolution available from an earth-rbased 
telescope. On the left is a Lunar Orbiter photograph of the same area 
which has been selected to show the resolution available from the same 
size telescope outside of the atmosphere. Further, elimination of the 
background glow associated with the atmosphere would permit longer 
exposure times in celestial photography, effectively enabling astronomers 
to see deeper into space. Realization of these advantages is contingent 
upon solving several technological problems. 

One of the more significant of these problems and the one which 
motivated the present investigation is the fabrication and maintenance 
of the primary mirror surface to the tolerance required for diffraction- 
limited performance. For the purpose of this paper, diffraction-limited 
performance is. defined as being obtained when the rms figure error 
over the mirror surface is less than a fiftieth of a wavelength, which 



k 



Figure 2.- Two views of the lunar surface indicating the increased 
resolution available from an orbiting telescope. 
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at 6328 A is l/2-|_i inch (l.27 x 10“^ cm). It is extremely difficult 
to achieve this accuracy with the monolithic mirror normally used in 
telescope applications as uncorrected thermal gradients, fabrication in 
a one "g" and operation in a zero "g" environment, and material reaction 
to stresses introduced during figuring all tend to cause distortion of 
the mirror surface. 

An interesting approach' which has been investigated as a means of 
solving the problem involves measuring the deviation of the mirror 
figure from the desired shape, generating the necessary control signals, 
and applying these signals to physically align the mirror to the desired 
shape. This concept, which is shown schematically in Fig. 3, has been 
investigated with both a segmented and a thin deformable mirror. The 
segmented mirror consists of a number of individual pieces or segments. 
This approach was selected because many of the effects causing surface 
deformations are reduced on mirrors of small size. In a Cassegrain 
telescope the incoming light is reflected by the primary and secondary 
mirror to focus behind the primary. In this case, the primary is also 
illuminated by monochromatic light which is returned to the mirror 
figure error sensor. The error signal is then processed and applied 
to the actuators which correct each segment in tilt and focus. 

This concept has been successfully applied to a small segmented 
mirror. ^ . However, the construction of large segmented primary 
mirrors from individually fabricated off-axis portions of a paraboloid 
(matched in focal length) would require new techniques of fabricating 
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diffraction-limited aspherics. In addition, the segmented mirror 
consists of a number of monolithic mirrors each subject to the same 
limitations as a monolithic mirror. . While these limitations are 
reduced in the relatively small segments, they are still of sufficient 
magnitude to require that the segments be fabricated from a substance 
which exhibits a very high degree of material stability. This is 
because the segments are corrected in tilt and focus only and any 
warping of a segment cannot be completely corrected. 

An alternate approach involves a thin, continuous -surface, 
deformable mirror that can be stressed into the desired shape by means 
of a large number of actuators arrayed across its rear surface. The 
greater control flexibility inherent in this approach shows promise 
of relaxing the material stability and fabrication tolerance require- 
ments. A laboratory model of a thin mirror is shown in Fig. 4 alongside 
a conventional monolithic mirror. The thin mirror is 30 inches ( 76.2 cm) 
in diameter and 0.5 inches (1.27 cm) thick. Actuators were located 
every 3-3 A inches (9*^ cm) over the rear surface of the mirror and 
were used to apply a corrective force distribution.^^ Preliminary 
operation of the system indicates that the actuators were able to 
reduce the initial figuring error to acceptable levels. This is 
illustrated in Fig. 5 which displays photographs of the diffraction 
pattern of the mirror before and during control. 

While this application appears to demonstrate the ability of the 
thin mirror approach to enhance the telescope performance, a n umb er of 
areas exist in which an improved design technique would be of significant 
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Figure 4.- Comparison of a thin deformable mirror with a mirror of 
normal thickness-to-diameter ratio. 
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Figure 5»- Comparison of the diffraction pattern of the uncontrolled 
mirror with that of the mirror during control. 




10 

value. In present design efforts the dynamics of the plant have teen 
largely ignored. This is a significant factor in the control of larger 
mirrors which have very lightly damped low frequency resonances. 
Selection and placement of the actuators is presently done on an 
arbitrary basis as a result of the present limited ability to relate 
the effects of these design decisions to system performance. In 
addition, the ability to more completely incorporate information on 
the disturbance characteristics, to the extent that it becomes 
available, is desirable. 

The purpose of this paper is to present a. general theory for the 
discrete control of a distributed-parameter system and extend the modal 
expansion technique to completely specify system performance. 

In the past little research effort has been devoted to the problem of 
obtaining discrete control of distributed-parameter systems and the 
results which have been obtained by Gould and Murry- Lasso^^ are limited 
to plants which have finite modal content. In this thesis the entire 
modal structure of the plant is considered for the problem of obtaining 
a specified level of performance while minimizing the number of control 
inputs to the plant. The method is demonstrated with examples and 
results are presented for a plant representative of a thin deformable 


mirror. 



II. MODAL EXPANSION TECHNIQUE 


A schematic representation of the plant under consideration is 
shown in Fig. 6 where L is a linear, time -invariant, separable operator. 
A typical example of a distributed system is the thin rectangular plate 
of Fig. 7 , where the deflection normal to the midplane of the plate 
w(x,y,t) results from the application of a transverse load density 
p(x,y,t). The equation of motion of the plate is given by'^ 


V 2 S x , y , t ) + p w(x,y,t) = p(x,y,t) 

&t 2 



where V 2 is the Laplacian in Cartesian coordinates, p the mass per unit 
area, and S = Yh.3/l2(l - v 2 ) is the flexural stiffness modulus of the 
plate. The deflection, w(x,y,t), is assumed separable, i.e., 

w(x,y,t) = c i (t)u i (x,y) . ' (V) 

Substituting 0=] into the homogeneous form of H yields, for a 
homogeneous plate of uniform thickness. 



+ flft L c 1 (t) = 0 



and 



i u i(x,y) = 0 
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p 

where is a constant whose value is specified by the boundary 
conditions. For the simply supported plate, which has boundary 
conditions^ 0 ) 


u i(o,y) = u ± (a,y) = u^x^O) = u ± (x,b) = 0, 


M 


and 


U. + V u. 

ax 2 1 ay 2 1 


= 0, 


x = 0,a 


OJ 


s X 3 

v — Ti u n + — r u n . 

ax 2 ay 2 


= 0, 


y = 0 ,b 


,(H) 


there exists' ' a denumerably infinite sequence 







where m, n = 1, 2, ... are the mode indices corresponding to the 
e igenf unction s 


u i = «nn(^y) = (sin m £ x)(sin n £ y) 


[5b] 


which form a complete orthonormal set (Appendix A). Consequently, 
the general solution to Eq.. M is 
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w(x,y,t) c i (t)u i (x,y) = ^ C 6 3 


i=l 


m,n=l 


and since the tu form a complete set the loading may he expanded in a 
uniformly convergent series of the form 

00 00 

p(x,y,t) -Ya^K^y) = ^ a mn (t)u mn (x,y) [6b] 


i=l 


m,n=l 


where 


li(-fc) =J J P(x,y,t)u i (x,y)dx dy 


C 6 5 


and T is the spatial region in which the plant is defined. 

Substituting [6a~j, [6b] , and OH into H and taking the Laplace 
transform with respect to time yields 

00 

^ (s 2 + o)^)c i (s)u i (x,y) = -^ a i (s)u i (x,Y) [j] 


i=l 


i=l 


where s is the Laplacian operator and, for convenience, the same symbol 
is used to denote a variable and its Laplace transform. Since the 
u ± (x,y) are independent the coefficients may be equated yielding 


Ci (s) = A i (s)a i (s) 


M 


\(s) = 


1/p 




S^ + 0^ 
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Based on Eqs. J^6a-cJ and jj3a-bj, the plant shown in Fig- 6 may be 
redrawn as shown in Fig. 8 which is a modal expansion of the distributed 
plant . 



Analyzer Synthesizer 

Figure 8-- Modal representation of the plant- 


More generally , Fig. 6 represents the functional equation 

w(z,t) = Lp(z,t) [V] 

where z represents a general spatial variable (in one or more 
dimensions) and the operator L operates on functions of time and 
distance. Laplace transforming Eq. {”$0 with respect to time yields 

w(z,s) = Lp(z,s) Q.oJ 

where the eigenvalues of L satisfy 


Lu ± (z) =A i (s)u i (z). 



Since the u^(z) are assumed to form a complete set, the separable 
functions w(z/s) and p(z,s) can be expanded as 







1 6 


w(z,s) = ^T C^sju-^z) 


i=l 


and 


CO 

p(z,b) =^T a i (s)u i (z). 


i=l 


Substituting OH- [l2a] , and [l2bj into O] yields 




00 00 00 

a i (s)u i (z) = ^T a i (s)A i (s)u i (z) (l3a) 
i=l i=l i=l 

0 

and consequently, 

•-(s) = A i (s)a i (s) [l3b] 

with a_^(s) defined by the transform of Eq. [6c]. Figure 8 is then 
the general modal representation of the class of distributed plants 
under investigation. In a function space where the eigenfunctions of 
L are used as the coordinates the system is represented by the infinite 
diagonal matrix 





In this reference frame, control of the plant output can be readily 
achieved by individually compensating each element of the diagonal 
matrix, as shown in Fig. 9. 
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Figure Decoupled compensation of A^(s). 


The significance of the modal control indicated in Fig. 9 is that 

it is relatively easy to express the system least squares performance 

in terms of the orthonormal modes of a vibrating structure. This is a 

phenomenon of particular interest in the mirror application since the 

integral square error is the desired performance index of an optical 
( 12 } 

surface.' ' With the error in the optical surface w e represented by 
the modal coefficients e^ the image index (expected integral square 
error) is given by 




e i(t) u i(z) 





p 

where E denotes the expectation and o\ is the variance of the error 

e i 

in the i th mode which is assumed to have a zero mean. The last steps 
result since u^z) is a member of an orthonormal set. Thus, the measure 
of image quality, J^, is a simple function of the variance of the mode 
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error which can he reduced by appropriate control one mode at a time 
as indicated in Fig. 9 • Relating the original signals in the system 
to those of the decoupled reference frame , the control system structure 
becomes that shown in Fig. 10 where e and a are column matrices whose 
elements are the modal coefficients e^ and , respectively. In 
practice the situation illustrated in Fig. 10 can only be approximated. 
The function of the analyzer is to determine the modal content of the 
optical surface error. The decoupled controller dynamics, represented 
by the matrix D(s) = diag C d i( s )] , is determined on the basis of 
standard design techniques (see Fig. 9 ) to achieve a satisfactory 
performance level. For a well-ground mirror the need for corrective 
action diminishes as the mode number, i, increases and control can be 
reasonably be restricted to the significant modes. The N controlled 
modes are denoted by the output and error N vectors c^ and e 1 ^ in Fig. 11. 
The finite (N X N) controller matrix is represented by D^. The 
function of the load synthesizer is to place an appropriate force 
distribution on the plate to correct for the modal errors in e^. Since 
the remaining modes are unmonitored (no corrective action taken), the 
ideal force distribution applied by the loading mechanism is 

N 

p ideal( z ^) =X a i (t) U i (z) C l6 ] 

i=l 

with none of the uncontrolled modes excited. Physically this corrective 
loading is applied by finite number of control manipulators which, in 
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the mirror problem, are comprised of displacement actuators in series 
with a spring acting against a backing structure. The spring is attached 
to the mirror by means of a pad intentionally introduced to restrict 
the excitation of uncontrolled modes (a point discussed in detail later 
in the paper). The actual (non-ideal) force density applied by these 
N actuators is given by 


N N 

p(z,t) Pj(z,t) aj(t)p^(z) [IT] 

j=l j=l 


where p • represents the force distribution resulting from the 
actuator, and the last step results under the assumption that each 
applied force distribution is separable in time and. distance. Expanding 
each of the Pj( z ) in terms of the eigenfunctions, Uj_(z), Eq.. [17] becomes 


DO 

p(z,t) =^T 


i=l 


N 

X 


j=l 


u ± (z) 


where 


[iSa] 


h io = / r h 

and comparison with M reveals 

N ' 

a i(t) 

j=l 
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In matrix form this relation between the mode force coefficients and 
the actuator signals becomes 


a = Ha 




where H is an ® x I matrix and is an N vector. Equation [19b] 
indicates the control elements excite all modes. Since only N of the 
more significant modes are controlled, Eq. [l9bj is partitioned as 



where is an N X N matrix, an «> x If matrix, a^ an N vector 

"D 

corresponding to the controlled modes, and a accounts for the remaining 

modal force coefficients. To provide the desired corrective vector 
N 

a the actuator inputs are given by 


a N , QjHj-l a N [ 21 ] 

where it is assumed that the actuator locations insure is non- 
singular - this point is discussed in detail in a subsequent section. 
Partitioning the matrix representing the plant dynamics into components 
corresponding to the controlled and uncontrolled modes, the overall 
system becomes that shown in Fig. 11 where disturbances q and q acting 
on the plant are included as equivalent displacements. When the mode 
number, i, is ordered with increasing frequency of vibration, = a^, 
the plant inherently performs modal filtering which attenuates the 
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higher modes so that their contribut ion . to the error (mirror distortion) 
rapidly becomes negligible as the mode number increases. Consequently, 
it is often assumed that only the first N modes are present, i.e., 

■D -p 

c and a are identically zero, and considerable simplification results. 
For example, let w^ be an N vector defined by the output at N different 
points. That is, 


w N = col w(z . ,t) 


O] 


where z. represents a measurement point. In terms of the mode 
displacement coefficients 


v 


N 


= U W c W 


where 


U» = 


u l( z l) u 2^ z l^ ' ’ * u u( z l) 

U X (Z 2 ) U 2 (Z 2 ) . . . UjjCZg) . ( 23 b] 


u l( z ])j) U 2^ Z N^ ’ * ’ U N^ Z N^ 

Under these conditions the mode analyzer becomes simply an operation 
on the N measurements; specifically 


c N = [/]\ N [24] 

where the sensors are located at positions to insure U W is nonsingular. 
The control structure of Fig. 11 then reduces to the N.'xH multivariable 



system^) shown in Fig. 12. This idealized representation, valid when 
the effects of the higher modes can he safely neglected, was derived 
hy Gould and Murray-Lasso and is treated in detail in reference 1. 
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Control of the low-order high-amplitude modes as indicated in Fig. 12, 

, ( 15 ) 

or in the decoupled form of Fig. 9 , presents the classic problem 
of controlling a resonant plant with a limited control effort (restricted 
actuator throw). While for large primary mirrors with low resonant 
frequencies this may he a substantial problem, in the present paper 
it is assumed that the disturbances q^(t) are slowly varying and of 
sufficiently small amplitude that any desired degree of control can' 
be realized. Negligible contribution' from the uncontrolled modes can 
generally be assured by permitting N to be arbitrarily large. However, 
in the present problem extremely accurate control of the optical surface 
of the thin deformable mirror for diffraction-limited performance is 
required with a minimum number of actuators. Under these conditions 
the effects of e^ and a^ are not negligible, but, in fact, represent 
the most significant system errors and the most important factors in 
evaluating design tradeoffs. 










III. ACTUATOR PAD SIZE AND LOCATION AS DESIGN FACTORS 


The uncontrolled modes enter the problem in two major ways. First, 
the actuators excite not only the controlled modes but, in general, 
all modes. To demonstrate, first separate the image quality index 
into two parts, i.e., 

where 



N 



i=l 



accounts for the error in the controlled modes and 


0 ° 

j r = Yj a h [ 25 0 

i=N+l 

accounts for the remnant error of uncontrolled modes. As larger 
actuator displacements are commanded in order to reduce to smaller 
and smaller values Jp. increases due to the effect of a^. Secondly, 
uniess direct measurements of the modes are made, a limitation on the 
ability of the displacement sensors to obtain an uncorrupted estimate 

■D 

of the N controlled modes results from the presence of e . 

Actuator Pad Size 

The function of the actuators, as indicated in the previous 

section, is to apply forces to reduce e^ while minimizing the 

excitation of the uncontrolled modes, i.e., ideally 

26 
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# nonsingular . 


00 


In the mirror problem the control manipulators are modeled by a 
displacement actuator working against a backing plate and a relatively 
soft spring which is attached to the mirror by means of a pad. The 
ability to approach the situation on Eq. GO is governed by pad size 
and location, which are factors under the influence of the designer. 

To illustrate the effect of pad size, consider the rectangular plate 
of Fig. 7 . The pads are assumed to be rectangular in shape and 
located as shown in Fig. 13* 



Figure 13*- Pad shape, size, and location. 
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The amount of force applied by the actuator is determined by the 
product of the spring constant K and displacement actuator position. 

This force loads the plate as indicated by Eq. [if] with 

Kf j(x,y) 

Pj(x,y) = 

f j(x,y)dx dy 




where f^(x,y) is the distribution of the force and a^(t) is the control 
input. Considering a force distribution that is constant over the area 
of the pad, Eq. [l8b] yields 



K 

(Ax) j(£y)j 



area 


mmx . nity 

sin sm dx 

a b 


hy 


. 


— sin 
ab 


mjtx. 

d 

a 


sin 



mit(Ax) ^ 

sin 

2a 

mty 

J 

b 

2a d 


nn(Ay) 

sin 


2b 


2b 0 



where i indexes the mode m,n. For this special case with constant pad 
size (i.e., (Ax) ± = (Ax)j = and (Ay) ± = (Ay) j = Eq. {[28]) may 
be rewritten in matrix form as 


H = GU’ 


[29a] 



29 


where U is the N x <» matrix, 


U = 


ui(zi) u 2 (z L ) 


u i( z l) 


U 1^ Z N^ U 2^ Z H^ ’ ’ ’ U i^ Z N^ 


[29b] 


the prime denotes the transpose of a matrix, denotes the point Xj_, 

yi> 


G = diag g i (A a ,A b ) , 


[ 29 ] 


and 


§i §mn ^ 


mnA 


n)t A 


sxn 


sm 


2 a 


2b 


mjtA 


2 a 


2 b 


O] 


Figure 14 contains a plot of (sin £ )/| . Assuming the controlled modes 

are m < and n 5: and pad dimensions are A a ^ ^ — and 

max 

. b 

, the maximum value of the argument for one of the controlled 

n max 

modes is {■ = it/2 which occurs when adjacent pads touch. The attenuation 
of the higher order modes by the pad is apparent from this diagram as 
the elements of G decrease rapidly for i > N and, in turn, decrease 
the output levels of approaching the idealized condition of Eq. [ 26 ], 
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Figure lk.- Plot of sin illustrating the filtering action of pad. 

The prefilter action of the pad is complemented by the transmission 

properties of the plate itself. For the rectangular plate of Fig. 6 

the relation between the applied loading and displacement output for 
"th. 

the i mode is 


c i(s) 


a i ( s )/p 

__ _ 

s + es + <x£ 



where a small amount of damping has been included. In response to a 
step input the steady transmission is 



O] 


lim 
t o° 
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where aS is given by Eq. O]. Assuming a 0.5 inch (1.27 cm) thick, 
30-inch (7^.2 cm) square plate with a Young’s modulus of 10 ^ pounds 
per square inch (70*3 X 10 gm/em^) and a Poisson's ratio of 0.2, 
this factor is 


lim _ 7.73 X 10~ 2 

t ->» a. / 2 , 2\2 

1 (m + n ) 



A byproduct of the pad’s desired effect on is a decrease in the 
transmission properties of H^. As a consequence, an increased effort 
is required to deflect the surface. This is readily demonstrated for 
the simply supported rectangular plate where 


e n = gV ■ (V0 

and with appropriate actuator placement is orthonormal (see 
Appendix B) . Consider the expected value of the norm of the actuator 
displacement vector given by 


N’ N 
a a 


= [a N '][G r ]- 1 [G H ] 


,N" 


1-1 N 
a 




Thus, as the elements of G- are decreased, the required control 
displacement and force increases. Since the plate itself was shown 
to perform substantial filtering, a compromise on the final pad shape 
and size is nomally employed. 
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Actuator Location 

Another critical factor on the H N and H R matrices is the placement 
of the actuators. The actuator locations must "be selected to prevent 
the occurrence of a singular II R matrix. Since the actuator locations 
most directly affect the matrix, knowledge of the mode shapes will 
generally permit the selection of the location of the actuators. This 
situation is detailed in the examples. Evaluation of the determinant 
may he used to verify the invertibility of I i . In addition, if a row 
of # is zero, the mode corresponding to that row cannot he excited. 

This is approximately achieved for most plants hy locating pads at the 
zeros (nodes) of the mode in question. For the rectangular plate, this 
is exactly achieved since the influence of pad location of H is 
delineated by U', see Eq. 03 - When the modes are ranked in order of 
importance, the desired goal is to control the first N and. null the 
next highest modes; however, this is not usually possible and trade- 
offs are required. For example. It may be possible to control a set 
of modes that are not the N most significant but be able to preclude 
excitation of the next highest modes or, alternatively, control the N 
most significant modes but not have the ability to preclude excitation 
of ' the most dominant uncontrolled modes. Since the controlled modes may 
be reduced to any desired level at the cost of some increase in the 
amplitudes of the unmonitored modes, the contribution of these higher 
modes to system error represents the most critical factor in system 
design. Determination of the trade-offs in actuator size and placement 


P 



is obtained through evaluation of the system performance index 
which is discussed below and illustrated later in two examples 



IV. STATIC PERFORMANCE 


The design objective is to minimize the image index defined as the 
expectation of the integral square surface error. As shown in Eq. O] 
this is a function of the mean square values of the modal error 
coefficients. Determination of Jj as given by Eqs. j^25a-cj is dependent 
on the nature, particularly the spectral content, of the disturbances. 
For the application to the control of a deformable primary mirror of 
an orbiting telescope, it is anticipated that the primary error sources 
will be initial figuring errors and relatively slowly time-varying 
thermal gradients. In this context it is reasonable to expect that tie 
system will generally be performing at or near its static values. 

In Fig. 11 the surface deformation due to the disturbances is 
defined in terms of its displacement modal expansion coefficients, q. 

No loss in generality results frctm considering the disturbances to be 
displacements since equivalent force distributions could be assumed. 

With reference to Fig. 11, and with ct^z, s) zero, the error in the 
control modes is given by 

e N = - (j + aV] q N . £ 56 ] 

For the static situation this reduces to 


S± " 1 + K, 


0 < i < N 


(XI 
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where is the loop gain (type 0 system) for the i^* 1 mode. When the 
loop transmission contains a pure integration (type 1 system), -» <»• 
Thus the error in the controlled modes can theoretically he reduced to' 
any arbitrarily small value. The expected value for J w then becomes 

[ 38 ] 



where cXq^ is the rms value of the static (or slowly varying) disturbance 
The' error in the uncontrolled modes is given by 


e E = - qK - tfj 


M 


= - q R + a V[hK]"V[i + aVIV 


[ 39 ] 


which under static conditions becomes 


N 




N + 1 < i 


j=l 


D*°5 


where \Jr. . is an element of the matrix 


+ = A E E E [a% K ] . 


M 


Assuming the modal coefficients of the disturbance are uncorrelated, i.e 


JO i f j 


\ 1 = J 
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the expected error in the -uncontrolled modes "becomes 


J R 


- I + l 


i=N+l 


j =N+1 


N 2 


- J R + J R 

o c 


00 


where J Rq is the value of the disturbance error in the higher modes 
without control. The second term Jr c is clearly positive and represents 
the increase in Jr that results from the actuator displacements 
required to control the errors in the first N modes. Since J Rc is 
finite, the series converges and the order of summation may be reversed, 
yielding 


where 



2 

9i = 


. 


[43b) 


' 3=N+1 


is a constant dependent on the design factors of actuator pad size and 
location as well as the natural modal filtering performed by the plant. 
Combining Eqs. jl?5a-c^] , [ 38 ^] , and 0*5] yields as the imaging index 

N r 


J I = J R 0 + 


i=l 


1 ''"a 2 

On 


1 + K* 


<Pj K i 

1 + K-i 


2 1 

\ 2 




The first term is unaffected by the control action, the second decreases 
as a result of the control action and the third increases. 
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Since Jr is constant , minimization of Eq. M with respect to 


Kj_ yields optimum performance when 


1 ■ 
K i 

cp? 




and 


N 


H 


L opt 


vI(rTK-j4 + J H 0 = Z 

1=1 i=l 


i=i 1 + * 


2 

<Pi 2 

a q n - + J R. 


2 y±- 


M 


with the controlled and uncontrolled. components given by 


N 


N 


J N 


opt 


i=l 


r2 _ 


2 

^i 




f-t \i + cp y ^ 

1=1 1 


) 


<£. 


O] 


and 


N N 2 

Jr -Y — cr 2 + J R S — cr 2 + J R . [kQ] 

R opt Z-j (i + k ) 2 R o Z_j (1 + cp2)2 ^i R o *“ 

i=l i 


i =1 (x + cp2): 


N 


Recalling that without control = > a n , it is seen that for large 

w o h_j ^i 
i=l 


loop gain the error in each term of JR 0 p t due to control is approximately 
l/K^ while each term in Jw 0 p^ is reduced by an additional l/K^. 

Thus, if significant improvement is to be gained in the optical surface 
by the above method, N must be selected sufficiently large and the 
actuator size and placement such that cp? = l/K. « 1. In the mirror 

■I* X 
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situation J p is negligible and Jl represents the major source of 
"o he 

concern. 

Because Kj_ is large, a type 1 system is normally employed and 



i=l 


The required value of K to achieve the desired rms figuring error can 
be minimized by the selection of pad size and location whose effect 
is manifest through the parameters cp?. 



V. MODE ESTIMATION ERROR 


The preceding analysis assumed ideal measurement of the controlled 
modal variables e^; however, in many applications it is neither 
practical nor possible to obtain direct measurements. In these cases 
an estimate is often derived from a spatial sampling of the distributed 
output. This is the case in the mirror problem where the most commonly 
used measurement of the optical surface is performed by the interference 
method illustrated in Fig. 15 - This mirror figure error sensor is 
a modified form of a Twyman-Green interferometer. In this interferometer 
two plane wavefront beams are formed from a common coherent source. 

One beam is reflected from a reference flat while the second is con- 
verted to a spherical wavefront whose center of curvature is that of 
the mirror. This wavefront is returned by the mirror and forms an 
interference pattern with the reference beam which is focused on the 
N discrete individual sensors. Periodic motion of the reference flat 
produces a sine wave of identical frequency at each detector. This 
converts the error determination from an amplitude to a phase measure- 
ment and permits the required sensitivity to be achieved. 

Under conditions where no modes except the first N exist, the 
relation of the modal coefficients to the N measured values is given 
by Eq. which for the mode error is 

e meas = OT 4* (V). 
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Figure 15 •- Schematic of the mirror figure error sensor. 




However, the presence of the high-order modes deteriorates this 
measurement since the actual measured vector is defined by 


w^ = U*e* + #e R 


where bP is the Ux ® matrix 


O] 


U 1 * = 


u N+l( z l) 


u- 


■N+j 


(zi) 


u N+l( z n) 


u n+o( z n) 


Consequently, the estimate of the modes in Eq. 0»] becomes 


-1 


-1 


e N = tW + rVH U*e R = e W + e N 

meas b. J v_- -J 


D*0 




The manner in which the measurement error e R , defined in Eq. [ 55 ], 
evolves is shown in Fig. l6. The disturbance error in the first N 
modes can be controlled to an arbitrarily small value (see Eq. Mi, 
while errors in the remaining modes cannot be counteracted. A major 
effect of the measurement error in Eq. is to introduce an 

additional error in the controlled modes. To illustrate, the vector 
e N is seen from Fig. 1 6 to be given by 


e N = - AW - q W 




oP = [ H N]-V[eN + el Q, 


[>h] 


where 










k3 


and substituting into O&a] and rearranging yields 

e N = - [l + aVJ\ N - [l + [ 55 ] 

For a type 0 system -under static conditions, Eq. Q?5~] becomes 


e - _ qi _ Ki e 
i 1 + K, 1 + K. i 


i — N 



or, for a type 1 system 

<!>» = - €l» = . [uH]‘VeE. [56b] 

Attention is now given to the error defined by Eq. [ 56 b] since, as 
previously established, most practical systems would possess an 
infinite loop gain. Two distinct cases are now considered: one where 

e R is dominated by the disturbances acting on the plant, i.e., 


JR 


R 



and, secondly, where the error in the higher-order modes contributed 
by the disturbance vector, q^, is negligible but the error introduced 
by the control effort is significant, . i.e. , 

e R = - A E H K a N = + A R #[a% N jj q N - 

= (x + c^VyV. [58] 

For ease in later calculations it is assumed that is negligible 

compared to I and that is nearly orthogonal permitting the contribution 
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to Jjj arising from to be determined from the norm of the measurement 
error at the N measurement points w^ = I^e* = U 1 ^. That is. 


J N 

meas 


. e[€ K cN] . b[v» 




where 



case I 
case II 



If the modal coefficients of the disturbances are uncorrelated then 
Eq. Q>9a] for case I becomes 


oo r N 




teas 




i=N+l n=l 


V 


[60] 


N 

V 2 

The sensor locations are chosen to minimize ) u i( z n) f° r as man y °f 


n=l 


the more significant high order modes as possible. This result is 
consistent with that concluded with regard to actuator placement and, 
consequently, the criteria for actuator and sensor placement are 
identical. Since increasing the number of sensors is relatively 
inexpensive, in many cases it will be desirable to have more sensors 
than actuators (or controlled modes) . If M > N sensors are used, a 
parallel development indicates that with increased measurements 



r N 




meas 


* 1 - l 


i=M+l 


n=l 
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yielding an expected improvement of 



[62] 


For the second case, of which the mirror problem is typical, the 
distortion in the higher modes is caused by the control action and the 
measurement error at the n^* 1 sample point is 



2 

From the properties of the disturbances, the expected value of w £ (z n ) 
can be determined and is given by 


N 


E[vf(z n )] = Y, I 


i=l 


u j( z n)*ji 


Li=N+i 




[64a] 


and 


i ' l meas 


N r N 

r 2 ( Z ) = y < y 

e v n y /A 

n=l i=l ^n=l 


w„ 


u .(z ) ijr . . 
3 n y ji 


j=N+l 


[64b] 


This error is on. the order of magnitude of the error in the uncontrolled 
modes Jp and represents the effect of estimating the modal coefficients 
from output measurements. Increasing the number of sensors, as suggested, 
will substantially reduce this error if the remnant is dominated by 
the first few higher modes. 



VI. SUMMARY OF DESIGN CONSIDERATIONS 


The performance index for the system is broken into two parts: 

(l) that contributed by the controlled modes, and (2) that due to the 
uncontrolled modes. Application of active control reduces the error in 
the controlled modes to' any arbitrary level while the disturbances 
producing errors in the (uncontrolled/unmonitored) higher modes cannot 
be counteracted. Further, the corrective forces applied by the finite 
number of discrete actuators excites additional errors in the higher 
inodes. When mode estimation is employed, errors due to measurement 
uncertainty are introduced into the first N modes. Thus, the total 
system error is given by 

= + X + Jr c 

where J w , Jp , and Jp are defined in Eqs. I 59a 1 and I 42 I. It was 
nneas n o n c 

observed that these errors may be minimized by factors under the 
designer's control. These factors are actuator location, pad size, 
and sensor location. Selection of actuator location permitted 
minimization of the excitation of the more significant uncontrolled 
modes. The pad size was seen to act as a filter which attenuates the 
effect of the control input in exciting the higher modes. The pad 
size is selected to cause the modal content of the applied force loading 
to drop off quickly above the II th mode. Combined with actuator locations 
that minimally excite the first few (most significant) higher modes, the 
pads together with the plant provide the desired modal filtering for the 
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remaining high-order modes. Finally, if estimation errors are to be 
minimized, sensor location requirements become the same as those for 
actuator placement. In some instances, additional sensors may be used 
to reduce the errors in estimating the modal coefficients. 

The above points are illustrated most clearly by the example of 
a simply supported, thin, square, flat plate with the following 
parameters: 

Thickness 0*5 inch (1.27 cm) 

Length 50 inches (76.2 cm) 

Width 50 inches (76.2 cm) 

Young's modulus 10^ psi (70-5 X 10^ gm/cnr) 

Poisson's ratio 0.2 

Table 1 contains values of oo? which are inversely proportional to 
mode transmission as given by the steady values of the A matrix. The 
modes are ranked, from most to least significant, in terms of decreasing 
(l/o^) and control of the low-order modes is desired. With the modes 
ordered with m on the ordinate and n on the abscissa, contours of equal 

n 

(oi^) become circles in the n,m plane. These circles may be approxi- 
mated by squares; for example, if 25 modes are to be controlled they 
would include 1 — n — 5 and 1 — m — 5 as shown in Fig. 17a. The next 
highest mode in each direction is m = 6 and n = 6 which has 25 mutual 
nodes equally spaced over the plate, see Fig. 17b. Actuator placement 
at these nodes (which can always be achieved since the m + 1 mode has 
m nodes) makes orthogonal and precludes excitation of any mode for 
which m or n equal 6. Alternatively, the first N most significant 
modes may be controlled as illustrated in Fig. 17c. This is achieved 
by the actuator placement of Fig. 17d. The unexcited modes in the 




Table 1.- Values of or in 10° radians / seconds z for a flat square plate 















































































(a) The m,n plane for the 
actuator placement of b . 



(c) The fi,n plane for the 
actuator placement of d. 


Figure 17.- The locations in the m,: 
and unmonitored modes of 
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(b) An actuator 
configuration utilizing 
twenty-five actuators. 
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(d) An alternate actuator 
configuration utilizing 
twenty-five actuators. 

plane of the controlled, unexcited, 
a flat rectangular plate’. 
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latter case are those with m = 8 or n = 8. In any case, final 
determination of the optimal trade-off requires detailed evaluation 
of the performance index. 

Since data regarding the disturbances to which a primary mirror 
surface is subjected are not presently available, a disturbance profile 
characterized by a modal force coefficient with a standard deviation of 




4 


'plate 

area 


386.4 



in pounds per square inch, for all i, was assumed. This profile yields 
an rms figure error of 


F = 




[ plate 
area 


«« 50 x 10" 6 inches 



for the uncontrolled surface which is in reasonable agreement with the 
figuring errors of the mirror whose diffraction pattern is given in 
Fig. 5* Pad size was selected to be 0.5 inches (1.27 cm) x 0.5 inches 
(1.27 cm). Fig. l8 displays the rms figuring error for a type 1 servo 
versus the number of actuators for placement of the type illustrated in 
Fig. 17b. The type 0 error was evaluated for optimal gain but did not 
provide a significant improvement for the error profile considered. 

The details of the procedure used to obtain Fig. l8 are contained in 


Appendix C. 
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The preceding sections contain the development and summary of 
design considerations for the discrete control of a distributed- 
parameter system. A simply supported flat rectangular plate has been 
used as an example because it possesses unique properties which clearly • 
reveal the results of design decisions which are obscured in most 
distributed system control problems. The following • section presents 
the application of the design technique to a plant which is representative 
of a thin deformable mirror and whose complexity is more nearly 
commensurate with that of plants generally encountered in practice. 



VII. MODAL CONTROL OF A FREE CIRCULAR PLATE 
Modal Representation 

The equation of motion of the free circular plate of Fig. 19 under 
forced vibration is 


2 2 . . <¥^w(r, 0jt) . . 

V S V w(r,0,t) + P 


[ 68 ] 


where w, p, and V 2 are expressed in cylindrical coordinates. Assuming 
solutions separable in r, Q, and t C •e., w(r,0yt) = c i (t)f i (.r)v i (e)j], 
the analysis follows that of the rectangular plate to yield 


where 


C ± (t) g 


dt 


2“ + £U L c i( t ) = 0 


[69a] 


a 2 v,(e) 


d 0 


2 + n ‘ \( 0 ) = 0 


[69b] 


r 2 d 2 Q , i (r)^J df i (r) 


+ r 


dr 


dr 


+ (+kfr 2 - n 2 )^) = 0 [69^ 




[ 7 °] 


Since the mode shape given by [69b] is periodic in 0, n is an integer and 
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u ± (r,0) = A[j n (kir) + BilnCkjrfJcos (n0 + 0 n ) (j 1 ] 


where A and 0^ are arbitrary constants, J n and I n are the n*'* 1 order 

Bessel and modified Bessel functions, respectively. The values of kj. 

and B^_ are determined through substitution of M into the boundary 

conditions, which for the free plate, arise from the absence at the 

free edge of both the bending moment in the radial direction and 

( 10 ) 

vertical shear j i.e.. 




r dr 


free edge 




and 

- o. [7a] 

free edge 

Under these conditions there exists a denumerably inf inite sequence 
of eigenvalues 

h -'*■> [ 73 ] 

for which the associated eigenfunctions form a complete orthogonal set 
permitting both p(r,0,t) and w(r,0,t) to be expanded in a uniformly 
convergent series as assumed in Eqs. [l2a] and [12b] . 

Table 2 contains values of k^ for several modes. Because of the 
importance of the mode shapes relative to design decisions, the radial 


( a I I 2 .1 a n 1 ] n 2 (v - 1) a.f c iW\ 




883 1.988 2.089 2.189 2.286 2.385 2.482 2.578 2.670 2.765 
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components of the first 21 modes, ranked in order of increasing 
frequency, co^, are plotted in Pig. 20. As a result of the nature of 
the 9 variation in Eq. 03 . each pair m,n is associated with two 
distinct modes given hy 


mode pair m,n = 


f mn( r ) cos ne 
f mn( r ) sin n0 


03 


for n f 0. For n = 0 a single distinct mode exists for each pair of 
m,n. To minimize the expectation of the square surface error, the 
actuators should affect control on the most significant modes as 
determined by the transmission factor 1 /pair and the disturbance profile. 
For the purposes of this paper, and as in the case of the rectangular 

p p 

plate, a force distribution with = 0 % is assumed such that an 

6 j, 

uncontrolled rms figure error of 50 x 10“ D inches ( 1.27 x 10“ cm) 
results. With this white modal disturbance the system objective is to 
exert control on those modes with the smaller values of k^n in Table 2. 
Actuator Size and Placement 

The actuators are again modeled as displacement actuators in 
series with a spring which is attached to the plate by means of a pad. 
The springs are relatively soft to make the effects of mirror 
displacement feedback negligible as discussed in Appendix D. The pad ' 
shape is a portion of an annulus bounded by constant increments in 
radius and angle. The elements hj^ are evaluated as 
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Ar^ 


h ij 


r j + 


0 th pad 


pad 


I II 


area area 


r j ‘ 2 


A9 i 
0 j + — l 


A Q . 

0, i 

J 2 


C J n( k an r ) 


+ B.I (k r)~l cos (n0 + 0 )r dr d0 
x n v mn 'J ' n' 


&Q 


for which the 0 dependent component is 


A0 


sin (nA0j/2) 


— t— — cos (n0. + 9 ) 

0 (nA 0 -i/2) J n 


where normally the increment A0 is constant over all j permitting the. 
inclusion of this component in the decoupled plant dynamics.- The r 
component is somewhat less tractable. Under the substitution 




[76] 


the r dependent portion of hy becomes 


•^mnr 

i r x 


r 3 + T", 


^mn 


Ar 

‘j " ~ 




l ) 


A [j n (]) + BI n ( T] 0 T) dTJ 


which is integrable if n is an even integer, but requires numerical 
integration or use of tables if n is odd. The effect in either case 
is that the radial component of hy decreases as k rnn increases. 
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In determining the placement of the actuators note from Fig. 20 

that each mode m,n has nodes at m distinct locations along lines of 

constant 0 and 2n nodes circularly. The placing of 2n actuators 

circularly at equal intervals results in their position coinciding 

with the nodes of one of the modes described by Eq. OI Consequently , 

this mode is not excited while the one spatially shifted 90 ° is. 

To determine the desired actuator placement a knowledge of the 

disturbance profile is required. Under the earlier assumption of a 

white modal disturbance spectrum (of = ), the objective is to 

i j 

control the modes with the smallest values of k mri . Inspection of 
Table 2 indicates that lines of constant k mn tend to form triangles 
connecting m to n where 

n = 2m . [77] 

If the controlled modes are n S. n max and m S rr i^y the controlled area 
of m,n plane is a rectangle (see Fig. 21a) which should approximate 
the region of the first N significant nodes. Actuator placement would 
fall at the nodes of the next highest modes which require 

^actuator = 2 ( n max + ^ Wax + ^ CW) 

actuators corresponding to the mutual ( m rnfly + l) nodes radially and 
2( n max + circularly. The controlled modes are the (2n max + l) 

Wax + bounded b y the rectangle m^, n^y plus the + 1 

excited modes for which n = n^y + 1. Thus, the total number of . 
controlled modes is H - 2(0^ + 1) (m^ + 1) = N actuator . The modes 


f 
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not excited by this actuator placement are those at n = ^max + 
whose nodes fall on the lines of constant 0 where the actuators are 
placed (those at n = n mflY + 1 spatially shifted, 90° are the extra 
m^y + 1 modes included in the N controlled modes). Additionally, 
the mode whose radial nodes are selected as actuator locations is not 
excited. The pertinent controlled and unexcited regions of the m,n 
plane are illustrated in Fig. 21a along with the corresponding 
actuator placement in Fig. 21b. However, because of the tendency 
of the lines of constant to form triangles as indicated in Eq. [77], 
control of an area in the m,n plane as indicated in Fig. 21c is 
generally desired. This can be accomplished by the actuator placement 
shown in Fig. 21d. Note that in the latter control scheme while the 
N most significant modes are controlled, it is not possible to preclude 
the excitation of the next most significant modes. The trade-off must 
be made on the basis of actuator location effect on the imaging index 

J I* 

While, at the present time, the imaging index has been detailed 
only for the rectangular control scheme illustrated in Fig. 21a, 
preliminary results indicate that control of the N most significant 
modes (e.g., see Fig. 21c) is preferred. The rms figuring error 
F = ^ J^/plate area based on the rectangular control scheme of Fig. 21a 
is plotted in Fig. 22 versus number of actuators for a plate with the 


following data: 
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(a) The m,n plane for (b) An ■actuator 

the actuator placement of b. configuration utilizing 

twenty-four actuators. 



(c) The m,n plane for the 
actuator placement of d. 



(d) An alternate actuator 
configuration utilizing 
twenty-four actuators. 


Figure 21.- The location in the m,n plane of the controlled, unexcited, 
and unmonitored modes of a free circular plate. , 
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Thickness 0*5 inch (1.27 cm) 

Diameter 30 inches (7 6.2. cm) 

Young's modulus 10' psi (70-3 X 10' gm/cm 2 ) 

Poisson's ratio 0.2 

These results indicate satisfactory control, yielding diffraction- 

limited performance, as defined by an rms figure error of less than 

0.5b inches, is achieved with less than 2k actuators. This is 

significantly less than the 6l actuators used in the present laboratory 
(lk) 

model which was determined by sectioning the mirror into equilateral 
triangles 3-75 inches on a side (the 3*75" x 3*75" x 3«75" x 0.5" thick 
triangles represent a thickness-to-area ratio near that normally found 
in monolithic telescope mirrors) . 



VIII. CONCLUDING REMARKS 


The modal expansion technique has "been applied to the problem of 
correcting and maintaining, to the tolerance required for diffraction - 
limited performance, the optical figure of a plant representative of 
the primary mirror of an orbiting astronomical observatory. The modal 
technique has been shown to be particularly appropriate for this problem 
by virtue of its relevance to a useful measure of image quality, its 
ability to decouple the system dynamics permitting simple control 
techniques to be applied, and by the extent of the insight the technique 
affords into engineering design decisions. 

For distributed plants subject to extremely accurate control, it 
is necessary to consider the effects on system performance of all of 
the modes - not just those which are subject to control. In fact, with 
the error in the modes under control reduced to any desired level, the 
major system error was shown to reside in the uncontrolled higher- 
order modes and this is increased by the control effort applied to the 
lower modes. For this reason the most significant design decisions 
are those related to the effects of the corrective control forces on 
the higher-order modes. The analysis presented in this paper describes 
the effect of actuator size and location on system performance , factors 
most critical to efficient design. The requisite conditions for 
minimizing the number of discrete control inputs required to achieve 
satisfactory performance were outlined and then illustrated in two 
design examples. The results for both the rectangular plate and the 
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free circular plate indicate that the thin deformable mirror can provide 
diffraction-limited performance; further, that this performance can he 
achieved with considerably less actuators than that required for a 
segmented mirror where the thickness-to-area ratio for each segment 
approaches that normally used in monolithic mirrors. 

The disturbance profile (if data on the effects of thermal 
gradients, spontaneous release of material stresses, or other factors 
producing distortion of the optical surface become available) can be 
readily incorporated into the design procedure. This is achieved by 
using the profile along with the transmission properties of the plant 
(plate) to determine the modal errors and the N modes yielding the 
largest errors controlled. Extension to more complex plants (e.g., 
shells), while requiring considerable computing effort, is direct. 
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APPENDIX A 


Determination of the Eigenfunctions and Eigenvalues of 
a Simply Supported Flat Rectangular Plate 

The equation of motion of a uniform plate in forced vibration^) is 


c^w(x,y,t) 

V^w(x,y,t) + p o — ■ - = p( x >y>t) 


St £ 




consider first the homogeneous equation. The modes of free vibration 
will be determined through the separation of variables technique by 
assuming 


w<x,y,t) - w 1 (t)w 2 (x,y) 
Sub stitut ing O] into [Al] 


[*€ 


+ p w 2 (x,y) 6 ''f ^ = 0 C A 3 

Ot^ 


dividing both sides by P 0 w p(t)w2(x,y) yields 

s v\f 2 (x,y) ± 


^l w l("t) 


P w 2 (x,y) W;L (t) dt 2 




Since each side of is a function of different variables, both 

2 

sides are set equal to oo^ - a positive constant. This yields the 
following equations 



Rearranging [A 5 ] yields 


[AT] 


(v 2 + a i\ffX^ ' “ > s|¥')" 2(x ’ y) ‘ 0 C a8 3 

The solution to[A8]is the sum of the solutions to each of the products 
Of or therefore to 

C«3 

[A10j 

is the Laplacian in Cartesian coordinates. Equations |_A-9j and 
CH are therefore 


V 2 + o^vJ^jwgCx^y) = 


^ " *\ltr 2ix ' J) = 0 


P 


r— 


■ <» 2 p w (x,y) 

( X ,y) - -i 2 = 0 

S 


—5 w 2 (x,y) 
dx 2 



= 0 





A further assumption is made that 


w 2 (x,y) = ^(xjw^y) 


[A12] 


In this case 


[All] 


becomes 


w^(y) 


d2 w^( x ) 


dx c 


+ W^(x) 


d 2 w^(y) 


dy 


+ CO. 

.2 “ 1 


\l 


w 5 (x)w4(y) = 0 (A13] 


Dividing by w^(x)w^(y) and rearranging yields first 


x a?, 3 (x) < x dg »hy) rr _ c 

w 5 (x) ^2 w^(y) d y2 ” Ms 


[A14] 


then 


w 


^(x) 


a %« ^ i p 

± CD.\ — = 


a 2 w^(y) 


dx £ 


iM S w 


4(y) 


dy £ 


[A!5] 


Since the left side is equal to a function of x, and the right a 
function of y both sides must be equal to a constant + p 2 . (The 

p p 

choice of sign on p is arbitrary since choosing as the constant - p 
will yield the same answers.) Equation r«ii becomes 


d^^x) 


dx 


+ 




= 0 
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The solutions to these equations are 


w^(x) = A^sin 


I - d 2 ± 




x + 



and 


w 4 (y) =A 2 sin (py + i 2 ). 




respectively. 

From equations QhlfQ, (~Al8j, and |”a 19] the solution to equation |a 5J 
is 


w 2 ( x ^y) = a 3 sin ( ^ 


-t. 2 ± 




x + sin (dy + £ 2 ) 



where A ^ is the product of Aq and A 2 - In order to evaluate the constants 
in jA2cfj the boundary conditions must be specified. For a simply 
supported flat plate they are^®^ 

w 2 (x^y) = 0 for x = 0,a [A2l] 

y = o,b 


% = 0 = - s 


k 2 v 2 ( x > y) s2w 2( x ^y) 


Sx 2 


+ v 


dy 2 


C A22 1 


x=0 

x=a 


and 
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/d 2 w 2 (x,y) 


S w„(x,y) 




Sy 2 


+ v 




y=o 

y=b 


[A23] 


In light of (j\2]~J the equations Q\.22j and |~A. 23 j may he written as 

= 0 

x=0 
x=a 

and 



S%(x,y) 

8x 2 


5 w 2 (x,y)| 

5y 2 


= 0 


y=0 

|y=b 


[A25] 


Using the conditions for x = 0, y = 0 from equation jA2lJ in [a2cT] 
yields 


fi l = *2 = ° 


[A26] 


Using the condition of equation [A2l] for y = b yields 


sin |ib = 0 


[A27) 


and therefore 


nit ~ n o 

M- - ">T n - 0^ 1^ 2 } . • • 


[A28] 



From the condition for x = a 


sin | - |f t ^ 


■li) 


a = 0 


[A29] 


or 


+ CD. 

^ “ i \l S 



p \ _ mjc 


= 20 - m = 0, 1, 2 


[a 3 o] 


Substituting jj\2&J and Qa3cTJ into Qv2cTj gives in [~A3l] the shape of the 
modes of free vibration, or eigenfunctions, of the plate 

w 2 (x,y) =A^sin S i n n,m = l, 2, 3, ... [A3l] 

The resulting mode shape is identically zero .for either m or n zero, 
consequently, equation [A31] is valid for the range of m and n 
indicated. 

That equation Q\3l] satisfies |a. 24J and [A25^j may be verified through 
substitution. Further, Aj may be chosen to satisfy the requirement 


p\> pa. 

/ / o)?(x,y)dx dy = 1 

“A UA 1 


0 u 0 


C«€ 


which yields 
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/ \ 4 mjoc . ruty 

w 0 (x,y) = -T- sin sm —■ 

2 X ' ab a b 

r -i s ^ wg(x,y) 

Since from [/tH = — Wg(x '^- 


s 

a 1 *- a 2 a 2 a 4 

W2(x,y) 

p 

a? + i? ay 2 + a y 2 

■w aa 

w 2( x >y) 






S 1 


2 2 2 2 


f4 4 _ _ _ . 

m « / x m jt n Jt , , 

-77 t + ” - — 5“ ^( x >y) 

”'La 


P w 2 (x 
+ “Vi w 2 (x,y)j 


CO? =S^ 

1 o 


2 .2 
a b 


2 2 
mr , n^ 

2.2 

a b 


= <4i 


This could, also have been obtained from (JI 27 J, jjV28j, and jll3oj. 
Since the solution to (6) is 


w 


00 

(x,y,t) = ) A cos(co + g ) sin 
v ' /_* mn N mn mn' 


nnx mjty 

sin — - — 


m,n=l 


A 36 " 


[a5?3 


»!<*) = *4 cosAjt + [a5S] 

The most general solution to the homogeneous form of M is 


[A59] 


Thus the eigenfunctions or modes of free vibration and their 
associated eigenvalues have been determined for a simply supported, 
flat, rectangular plate and are given by equations [A34^ and {A37^] 


respectively. 



APPENDIX B 


Determination of a Set of Actuator Locations for Which For 
a Simply Supported Flat Rectangular Plate is Orthogonal 

The purpose of this appendix is to indicate, for a simply supported, 

flat, rectangular plate a set of actuator locations which make the 

matrix (eq. 09) orthogonal. The equation of motion of a beam^ is 


h o 

2 5 w(x,t) _ S w(x,t) 


6x 


dt‘ 


[bi] 


The method of separation of variables is used, consequently, 

w(x,t) = w 1 (x)w 2 (t) [W] 

Substituting (? 2 ] into [Bl] gives, after dividing both sides by w(x,t). 


2 d^w^x) 1 d 2 w 2 (t) 


w- 


r l(x) W 2 (t) 


dt c 


[B3] 


Both sides are set equal to a constant 


d V( x ) k 


dx' 


7 wj_(x) = 0 


^ w 2(t) > 

+ 7 w 2 (t) = 0 


dt c 


(^B^a) 


Qb4t0 


The solution to [bW] is 


w 2 (t) = C 1 sin 7 2 t + C 2 cos ? 2 t 


I B 0 
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The solution to the remaining equation, QAa), will yield the modes of 
free vibration of the team. Since the exact solution will be used in 
the sequel this information will be obtained first. 

Equation is factored as 



A( c - 


/\dx £ 


+ 7 —hiW = 



and the solution is the sum of the solutions to 


[B7] 


and 


kdx 


4 - vW*) = ° 



= 0 



Thus 


w-,(x) = C^cos — x + Ck sin — — + Cc;cosh — — +C^sinh — — pB' 

9 \[r \F \Tr \pr L 

at this point the boundary conditions are brought in. For a simply 
supported beam 


w L (0) = 0 

[B1O0 

o 

ii 

47 
% — ' 
c— l 

IS 

(BlOb) 
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and there are no moments at either x = 0 or x = b 



= 0 

x=0 

x=b 



These conditions require 


0 = Cj + 



0 = C'zcos — — b + Ch sin — — b + Cc cosh —Lr b + Cgsinh —= b (pilb] 

; ^fr ' 4 \J7 ' 

0 = - ^c 5+ £ C 5 [Bile] 


^1 

1 

II 

0 

Gx cos 

7 

"b + Ck sin 

— b 

y2 

+ — 

CRCOSh 

b + CA sinh — — b 



\f¥ 

\F 

\[7 


\[r \[7 


[BUd] 

From equations 0511aj and 0Bllcj it is determined that 


c 5 = c 5 = 0 

leaving 



Vb y 

0=04 sin + Cg sinh — b 

\Tt \J~r 

[bu] 

7b 7 

0 = - C 4 sin — = + Cc sinh — = b 

\F 

[Big 


A nontrivial result requires the determinant of the coefficients of 
C4 and Cg in equations (?i3 and [b^ to be zero. Consequently, 
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^ • 7 b . , 7 b 

2 sin — — sxnh — — = 0 

\F \p 


[bis] 


or 


7 b 

n/T 


= mjt 


[bi.6 ] 


Substitution of [bi 6] into C b1 € indicates that Cg = 0 and 


7 _ mit 

\fr b 


[bit] 


The eigenfunction of the homogeneous simply supported beam is 


w^(x) =04 sin ~ x 


[bib] 


. The solution to the problem is now considered through the method 
of finite differences. A number of stations are located at equal 
intervals, l, along the beam as shown in the following sketch 


L 

]<- 1 ->l ^ 

, 

I 

□ 

0 

/ 

X 1 x 2 . • 

/ 

• • x n-l 

x n 

b 


and the equation 


dM 

"a? 


wN = 0 



is written at each point, where w N = col co(x 1 ),... The fourth 

derivative approximation used can be obtained by first obtaining the 




8o 


Taylor series expansions about x r of w (x r+1 ), ^x^), w ( x r _i); 
w ( x r-2)* These expressions are 

w ( x r +i) = w ( x r ) + w'(x r )Z + w "( x r)~ + • • • + w (n) (x r )~ [B20a) 

2 

w ( x r-l) = w ( x r ) - w '( x r )^ + w"(x r )- . . . + (-i) n w ^ n ^( x )— 

2! n! 

[b2O10 

w(x r+2 ) = w(x r ) + w’(x r )2Z + w "(x r )~p- + . . . + w^ n \x r ) ~y— 

[ B 2O0 

w ( x r-2^ = w ( x r^ " w'CXr) 2 ^ +'w"(x r )- | - p + . . . + (-l) n w( n ) (x r ) 

(B20d) 


subtracting four times (^B20aJ + |j320bJ ) from 0B2Ocj + |^B20dj yields 


4 w ( x r+l^ + w ( x r „]_) + w ( x r+2 ) = " 6w ( x r ) + 2 ~' T^~ 

clx 


.^(2 4 - "0 


X^QCj. 


+ 


+ 2 d2 M x ) 


dx* 


,2m 


2fl (a 2 ” -4) 
(2m)! 


x=x T 


+ . • • 


After rearranging equation [B2l] becomes 
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w (x r _ 2 ) - 4w ( x r _i) + 6w < x r ) " 4w ( x r+i) + w ( X r+2 ) d 4 w(x) 


dx 


x=x r 


z 2(m-2) 2 2(m-l) 

“ n . — “f* « « • 

(2m)! 

x=x r 

(B22] 

« 

The term on the left will he denoted L^. The finite difference method 
approximates the plant equation as 


+ 2 d2nV ( x ) 
dx 2111 


B fd w (*r> - %- ”< x r ) - 0 C B2 3] 

T 

Collecting the expressions for each point and arranging them as 
indicated by the definition of yields the matrix equation. 


kfdm 


X 1 

X s 

l- 1 

w 1 

2 

w( Xl ) 

• 

- 7fd 

• 

• 

• 

w (*n) 

T 

• 

• 

w(x n ) 


[B2*Q 


The solutions to this equation are the eigenvectors of the finite 
difference - matrix representation of the plant, Lf^. The matrix 
is a real symmetric matrix and consequently, has orthogonal 
e igenf unctions . 

If the right-hand side of equation [b 22^ is used instead of the 
left it is possible to obtain an analytic expression for the numerically 



82 


obtained finite difference answers. That is, consider the equation 


_ d w( x ) 

ax 1 * 


fd 


w(x) + . 


,2m , x ,2(m-2) o/ . 

+ 2 i_h*i ±1 : + . . . 

dx 2 " 1 ^ 2m ) ! 

[ B 2 5 ] 


The answer to this equation evaluated at x = x^, is equal to the finite 
difference answer. 

Consider the function u ^ = C4 sin — x which, from [B18] , is seen 
to be the exact eigenfunction. Substituting this expression into [B25] 
yields, after cancelling u^ 


7^ ^fd 

0 = — — + 

T 2 t 2 


. + + . . . [ b26 ] 

V7/ (2m).» L J 


Thus Uj_(x) is the solution to [b 20, and u^(xq) = w(x^), when 


fd 


= ^2 + 


\2m 


(-1)° / , y (2n 2 ( m - 1) , 

j2 


m > 3 


[B27] 


l is the separation between adjacent stations on the beam and it can 

be made as small as desired. Consequently, equation [B27] indicates 

More significant for the 
th 


7fd 7 

the resultant convergence of . to 

VT \pr 


purpose of this appendix is the fact that the i component of the 
j^k eigenvector of the finite difference solution is equivalent to 
the j'kk eigenfunction of the beam evaluated at a location corresponding 
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to the i^h point in the finite difference representation of the beam. 

Since the matrix of Eq. 03 is 


u i( x i) • • • u n( x i) 


u l( x Ij) • • • %( X K ) 


[B28] 


each column corresponds to one of the orthogonal eigenvectors of L ^ m 
and the matrix is itself orthogonal. 

To this point the proof has concerned the solution for a beam, 
while it is desired to show 


X u i( x k> u j( x k) = 0 1 I j [ B2 9D 

k=l 

for the plate. For the plate, Appendix A shows that 


, . m- mxir n • «y k 

U-! (x v ) = C-i sin — sin — — 

1 K 1 a b 


consequently, Eq. [B29] may be written 


0°3 


N N 

X u i (x k )u j (x k ) = X ^i^kKi^kKjK^jCyk) [B3l] 

k=l k=l 



jlt-ls 
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Consider an array of locations 


N = XY 


[B32] 


with X locations in the x direction and Y in the y direction. Thus 
becomes 


X Y 

a xi( x k) u yi^k) u xj( x k) u yj^ x k) = / / u xi( x A) u yi^yB) u xj( x A) u yj(YB) 

A 3=1 B4=l 

[B33j 

X Y 

°r = a) 

A 5 =l Bk=l 

Each component in the second expression for the right-hand side of 0 B »] 
is equivalent to the beam, consequently, either the first or second 
summation will be zero unless i = j . 

Equally spaced points will provide an orthogonal matrix for the 
purpose of relating performance specification in the original and 
transformed systems. 



APPENDIX C 


Evaluation of the RMS Figuring Error' for a Simply 
Supported Flat Rectangular Plate 

The purpose of this appendix is to describe in detail one of the 
procedures used in determining' the results contained in Figure 18. As 
indicated in equations O] and [st] the exact determination of Jj would 
require the evaluation of an infinite number of terms. In using a 
finite number of terms p to approximate Jj and, therefore the rms error, 
it is desired to select such a P, if possible, which would place a bound 
on the amount by which the approximate value of Jj would differ from 
the true value. 

The value of P which should be used is a function of the individual 
terms in the sequences 

a x , a 2 , ... a n ... [Cl] 

and 

«S> ••• ••• [02] 


Since 

a i 

Ci = -±££ 

SS 

one procedure might be, for monotonically decreasing values of Cj_ ss , 

to observe values of c-j_ ss in [cj] for increasing i until values of 

Cj_ are obtained which are significantly less than the accuracy desired, 
s s 
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The value of Jj may be determined for two values of P in this neighbor- 
hood to determine whether or not Jj has been obtained accurately enough. 

While such a procedure would be adequate the specific nature of 
the present problem permits the selection of P on a more rigorous 
basis. For a square plate the expression contained in Appendix A for 
the eigenvalues becomes 




From Eqs. Qfj, [15J, and and [c^j 


2 

a 8 a i 


'a ob v- “ 8 a i 

» (*,y,t)a* ay " + n aj4 




Since m and n each take on the values of all the positive integers 
the right-hand side of [c 5 ] may be rewritten 



00 00 



[i 2 + Ff 



In order to remove one of the infinite summations, use is made of the 
symmetry of the eigenvalues by writing C«Q as 
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where aj_j is considered a constant as has already been assumed (see 
Eq. [66]) . The terms i 2 + j 2 in the series 


oo i 

II 


i=l j=l 


1 

[l2 + 



are placed into correspondence with the positive integers in the order 
indicated by [c8]. It is desired to indicate for the i and j 
corresponding to the k^h value in [C8] that 




[cs] 


In the first n values of i there' are 

n 

T =I 1 

i=l 

terms where T is 

T = n t n + 

2 


[cio] 

[oil] 


2 2 

The j = 1 term is the minimum value of i + j for any value of i. 
Since, for i = n 


if 


n 2 + j 2 > n 2 


n 2 > n ( n + X ) 




2 
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inequality [C9j will "be obtained. Since [cig is valid for 


n > 1 


then ^C9_ holds and the right side of [ct] may he written 


0*0 


8 2 00 1 
a a U < 


Ji Ji fi 2 + 


8 2 ” 

< 2 V 1. 


O 5 * 4 ] 2 h k fr 2 + 


SjA k^ 

i=l 


[C15] 


The advantage of the last series is that its sum has a hound. The 
partial sum for the first 2^-1 terms is 


o 1 .1 .,1 . 1 

Qot i = — r + ~ r + t + • • • + r 

2^ 3 1 * (2 t - l) 1 * 


[C16] 


which is less term by term than 


Q 


2t_l 




which is 


thus 
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The second of the two terms is always negative and goes to zero as 
2 t - 1 increases. Thus the sum of the primed series from the 2^-1 
term to the end, Q , is less than the second term on the right side 

of [19] 




2 5 - 1 



Since the primed series is greater than by term than the k series, the 
sum of the k series over the same terms must also be bound by the same 
amount. 

, v t-1 


I 

k=2 t 


1 MV 
? 15TT 


[C21] 


The original series is smaller than the k series term by term, 
therefore, since 

2 —1 Q p 00 p p 

n b r~> a af v -1 a af 

- u i=l C Sir4 ] K + n i] 4 i= 2 t [s^] 2 [nf] 


[C22] 


then 


^a ob 


o w o 


2-1 g g 

w 2 (x,y,t)ax dy < Y -klk 


8 a il 


. 2/1 


t-1 


a 


23) 


i=l [s^] 2 [?i + n f [s^] 2 25-1 
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and 


n b r-i 

cd (x,y,t)dx dy - ^ 


2 t -l 


a 8 a? 

i 


< & 


2 1 
a i 3 
8 2^ 


t-1 


kt^y&L + »d k M 2 23 - 1 


O^] 


The inequality may be used either to place a bound on the error 
for a fixed number of modes considered in the evaluation, or conversely, 
it may be used to select the number of modes which must be considered 
to keep the error below a certain level. It should be noted that the 
value of a^ in 0*0 is the result after control and must be chosen 
conservatively enough to reflect the induced error as well as that 
already present. Once the bound on the error has been selected t may 
be determined from C«*D- This equation was based on the use of 2 ^ - 1 
terms in the series, consequently a square n x n array of modes cannot 
be used for which 

n(n + 1) > 2 t _ x 
2 

Finally, the tightness of the bound is dependent on how closely 
the k series approximates the true series for the terms after the 2^-1 
term. The bound may be tightened by observing actual deviations and 
adjusting the inequality [C 24 ] by the appropriate amount. 



APPENDIX D 


Mirror Displacement Feedback 

The actuators considered in this appendix are modeled as types 
which are being considered for actual usage^”®^. This model consists 
of a pure displacement actuator acting against a spring and a backing 
plate which is stiff relative to the spring. To obtain a specified 
force the displacement actuator is commanded to a new location relative 
to the undeformed mirror. 

If the mirror has deformed, the displacement of the mirror will alter 
the magnitude of the applied force. It is assumed, arbitrarily, that 
one form of this displacement feedback might take is 


Pi(x,y>t) = o^tJPjfoy) -Kw i (x,y,t) 0H] 

where p i (x,y,t) is the force density applied by the i* 3 * 1 actuator, and 
where 


W£(x,y,t) = w(x,y,t) 


over the i^* 1 pad area , and 


Wi(x,y,t) = 0 


elsewhere , and 


a-i(t) = 


KZ.(t) 


ff M 51 ' 


y)dx dy 




D»] 
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w^x^t) may be expanded as 


w j _(x,y,t) =2 j c ij( t ) u j( x ' y ^ 

j=l 



where 


c! .(t) 
i«T ' 



w i (x,y,t)u.(x,y)dx dy 



Since w^(x,y,t) is zero except over the pad area of the i^ pad C*Q 
can be written 



(t) 



w(x,y,t)uj(x,y)dx dy 


ith p a d 

area 



Substituting the modal expansion of w(x,y,t) yields 


c h (t) 


oo 

= J f Yj c k^ t ^ u k^ x,ar ^ u j^ x,y ^ dx dy [BY] 

. +v> , k— 1 


ith pad 

area 


Interchanging the order of summation and integration yields 


00 



k=l 


J J u k (x,y)uj(x,y)dx dy 


i^k pad 
area 
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The total feedback force, f s (x,y,t) is 


N 


f“ s (x,y,t) = - K ^ w i(x»y^) 
. i=l ' 


Substituting 0*3 into 03 


N oo 


f s (x,y,t) = - K ^ ^ cjj(t)Uj(x,y) 
i=l j=l 


03 


O 0 3 


Interchanging the order of summation 


f s (x,y,t) = - K ^ uj ^ c[j(t). 

j=l i=l 


[dh] 


By comparison with equation 03 f s (x,y,t) can be expressed modally as 


f s (x,y,t) 


OO 

-I 


a j (t)uj(x,y) 


[D12] 


0=1 


where 


N 

a j (t) * - K Z 

i=l 


03 


Equations [D13] and 03 can be used to put the expression for aj(t) 
into a different form 


_a i W 

K 


N oo 

= Yj x c k^^ r i u k( x ^) u o( x ^^ dy [p 1 ^ 

i=ife4 'iv, J . 


i"kk pad 
area 



Interchanging the order of summation 


00 N 

- 4 ^ -£*(*)£ J J %(x,y)uj(x,y)dx dy Q)15] 

k=l i=l ith pad 

area 


thus 


- a ' (t) = + KZ c(t) 



where Z is an » x « matrix which has individual elements of 


J j l 


N 

= Y_, f f uj(x,y)u 2 (x,y)dx 


i=l -,-th 


i 0ii pad 
area 


[D17] 


and where the spring constants are all assumed equal. 

The matrix Z is located in a local feedback loop around the 
diagonal plant matrix since it describes amplitudes in the force modes 
as a function of the displacement mode amplitudes. Since, in general, 
the elements of Z are non-zero, this feedback causes the system to 
become coupled. In this particular application the coupling exists but 
can be made to have negligible amplitude. Since the displacement of 
the mirror is on the order of microinches the spring constant is adjusted 
so that the actuator throw required is on the order of inches. 

Alternatively, a form of spring feedback may occur which can be 
treated without introducing coupling. Suppose that the expression for 
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Pi(x,y,t) 


K(ldj - Imjjpjfey ) 
J J^p’ 1 (x,y)<ix; dy 



which would he obtained by decreasing the pad area until the mirror 
displacement over the pad area is constant. This, depending on the 
manner in which the pad is bonded to the mirror, appears to be a 
reasonable assumption. Expanding 0Dl8^J yields 


Pj_(x,y,t) = 


KZ d,- Pi^y 1 * 


P i (x,y)dx dy 



KimpPil^y) 


P i (x,y)dx dy 


03 


This equation can be analyzed by a procedure similar to the preceding 
paragraph, or equivalent results may be determined from an inspection 
of the appropriate block diagram. Proceeding as previously, the second 
term in 03 is 


f s i ( x ^y^ t ) 


= - KVjPi^y) 



Pi(x,y)dx dy 



where ^ is the constant value of w(x,y,t) over the pad area. It is 
desired to express f s (x,y,t) in a modal expansion 


f s .(x,y,t) = ^ a!y(t)uj(x,y) 
j=l 




w.(x,y,t)B.(x,y)u (x,y)dx dy 

J [D22] 
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Since Wj_(x,y,t) is constant over the pad area and 3j(x,y) is zero 
elsewhere OD becomes 


a « (t) ' 


KZ m . (i) 


ith pad 


Pi(x,y)uj(x,y)dx dy [p2^j 


p i (x.,y)dx dy ar g a 


The integral has been previously evaluated as h (Eg.. [18b]) , therefore 


, x - KZ m .(t) 

a: .(t) = — hi • 

r r J 1 

/ / P-i (x,y)dx dy 


Qd24} 


and if the integral of B. (x,y) is equal to that of pi(x,y) then 


where 


aj(t) = 


- K 


HZ N (t) 


PiCx^yJdx dy 


[ D 25] 


J H (t) = 


w(xj,yj,t) 


[w(xj ; y w ,t)_ 


= w N (t) 


[D26] 


In this case the general result indicates that there is coupling in the 
system. If the system is assumed to have only finite eigenfunction 


content then 
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N 

wN = Yj c i ( t ) u i( x ^) 

i=l 

If [D2 7 ] is placed in matrix form 

w N = c(t) 

where 



[d 28 ] 


U J 


rN = 


Ui(xi,yi) . . . U N (x!,yi) 


Ui(x K ,y n ) . . . %(x N ,y n ) 


[D29] 


The spring feedback loop as described by [Dl8] and [D2 5 ] is shown in 
Figure D-l 



Figure D-l 

The spring feedback loop. 
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As indicated previously, under appropriate conditions, the matrix can 
be written as the product of a diagonal matrix and a nondiagonal matrix 
as in 3^ • In this case the spring feedback loop becomes that shown in 
Figure D-2 



Figure D-2 


The spring feedback loop for = ff'V* . 


Both G and 


are diagonal elements. Consequently, 


3 i (x,y)dx dy 


when u*’ = iT 1 the system is decoupled (see Fig. D-3) and the effects 


of the spring feedback can be readily included in the system design. 



Figure D-3 

The decoupled spring feedback loop. 











99 


In summary, under a specific set of assumptions, the effects of 
the mirror displacement feedback can be treated without introducing 
coupling effects into the system. In general, coupling effects are 
present, however, the mirror displacement feedback is rendered negligible 
through an appropriate choice of spring constant. 
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